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METAL ARTIFACT CORRECTION IN COMPUTED TOMOGRAPHY 



The following relates to the diagnostic imaging arts. It finds particular application 
in computed tomography imaging of a subject that includes high density regions such as metal 
implants, dental fillings, and the like, and will be described with particular reference thereto. 
However, it also finds application in other types of tomographic imaging such as single photon 
5 emission computed tomography (SPECT), positron emission tomography (PET), three-dimensional 
x-ray imaging, and the like. 

In computed tomography imaging, so-called "metal artifacts" arise when the 
imaged region of interest contains metal implants, dental fillings, or other regions of high radiation 
absorption. Typically, metal artifacts appear in the reconstructed image as streaks emanating from 
10 the high density region. Projection line integrals passing through the regions of high density are 
highly attenuated leading to substantial measurement errors, and the filtered backprojection or 
other reconstruction process translates these measurement errors into streaky image artifacts 
commonly known as metal artifacts. 

A previous method for correcting metal artifacts includes performing filtered 
1 5 backprojection to generate an uncorrected reconstructed image, identifying a region of high density 
in the uncorrected reconstructed image, and replacing projections that pass through the high density 
region with synthetic projection data having reduced absorption attenuation values. The corrected 
projection data again undergoes filtered backprojection to produce a corrected reconstructed image. 

This known method works well for certain imaging applications in which there is a 
2 0 single, well-defined high density region surrounded by much lower density tissue. It does not work 
well, however, with a plurality of high density regions, or where there are medium density regions 
in addition to the high density region. For such composite imaging subjects, metal artifacts are 
reduced but remain veiy visible in the corrected reconstructed image, especially between high 
density and medium density regions. In medical imaging applications, medium density regions 

2 5 typically correspond to bone while high density regions typically correspond to metal implants, 

dental fillings, operation clips (used in certain interventional computed tomography applications), 
prosthesis devices, and the like. Hence, in medical computed tomography imaging, the region of 
interest commonly contains medium density regions. 

The present invention contemplates an improved apparatus and method that 

3 0 overcomes the aforementioned limitations and others. 



WO 2005/008586 



PCT/IB2004/002325 



According to one aspect, a method is provided for producing a corrected 
reconstructed image from acquired tomographic projection data. Acquired projection data 
corresponding to a region are reconstructed into an uncorrected reconstructed image. Pixels of the 
uncorrected reconstructed image are classified at least into high density, medium density, and low 
density pixel classes. Pixels of the uncorrected reconstructed image that are of the high density and 
low density classes are replaced with pixel values of the low density pixel class to generate a 
synthetic image. The synthetic image is forward projected to generate synthetic projection data. 
Acquired projection data contributing to the pixels of the high density class are replaced with 
corresponding synthetic projection data to generate corrected projection data. The corrected 
projection data are reconstructed into a corrected reconstructed image. 

According to another aspect, an apparatus is disclosed for producing a corrected 
reconstructed image from acquired tomographic projection data. A reconstructing means is 
provided for reconstructing acquired projection data corresponding to a region into an uncorrected 
reconstructed image. A classifying means is provided for classifying pixels of the uncorrected 
reconstructed image at least into high density, medium density, and low density pixel classes. A 
pixel replacement means is provided for replacing pixels of the uncorrected reconstructed image 
that are of the high density and low density classes with pixel values of the low density pixel class 
to generate a synthetic image. A forward projecting means is provided for forward projecting the 
synthetic image to generate synthetic projection data. A projection replacement means is provided 
for replacing acquired projection data contributing to the pixels of the high density class with 
corresponding synthetic projection data to generate corrected projection data. The reconstructing 
means reconstructs the corrected projection data into a corrected reconstructed image. 

According to yet another aspect, A radiographic scanner is disclosed including a 
computed tomography scanner and a processor. The scanner includes at least an x-ray source, a 
detector array, and a rotating gantry on which the x-ray source and the detector array are mounted. 
The scanner acquires tomographic projection data during rotation of the rotating gantry. The 
processor produces a corrected reconstructed image from fee acquired tomographic projection data. 
The processor performing a method including: reconstructing acquired projection data 
corresponding to a region into an uncorrected reconstructed image; classifying pixels of the 
uncorrected reconstructed image at least into high density, medium density, and low density pixel 
classes; replacing pixels of the uncorrected reconstructed image that are of the high density and low 
density classes with pixel values of the low density pixel class to generate a synthetic image; 
forward projecting the synthetic image to generate synthetic projection data; replacing acquired 
projection data contributing to the pixels of the high density class with corresponding synthetic 
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projection data to generate corrected projection data; and reconstructing the corrected projection 
data into a corrected reconstructed image. 

One advantage resides in improved artifact correction for images having a plurality 
of discontinuous regions of metal or other high density material. 
5 Another advantage resides in improved artifact correction for images including a 

metal or other high density region and a bone or other medium density region. 

Numerous additional advantages and benefits will become apparent to those of 
ordinary skill in the art upon reading the following detailed description of the prefened 
embodiments. 



The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. The 
drawings are only for the purpose of illustrating preferred embodiments and are not to be construed 
15 as limiting the invention. 

FIGURE 1 diagrammatically shows a computed tomography imaging system 
including an artifact-correcting reconstruction processor. 

FIGURE 2 diagrammatically shows a preferred set of pixel density thresholds 
defining a preferred set of pixel density classes. 
20 FIGURE 3 shows a block diagram of the artifact-correcting reconstruction 

processor of the computed tomography imaging system of FIGURE 1 . 

FIGURE 4 diagrammatically shows sinogram representations of acquired, 
modified, and synthetic projection data at several points in the reconstruction processing. 

FIGURE 5 diagrammatically shows reconstructed and partially processed images 
25 at several points during the reconstruction processing. 



With reference to FIGURE 1, a helical conebeam computed tomography imaging 
scanner 10 includes a radiation source 12 that produces a radiation beam directed into an 
3 0 examination region 14. The radiation beam interacts with a region of interest of an imaging subject 
disposed in the examination region 14, producing spatially varying absorption of the radiation as it 
passes through the examination region. A radiation detector 16 detects the absorption-attenuated 
radiation after it passes through the examination region 14. 

In a preferred embodiment, the radiation source 12 produces a fan-beam or 
35 cone-beam of x-rays. The radiation source 12 and the detector 16 are preferably mounted in 
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oppositely facing fashion on a rotating gantry 20 so that the detector continuously receives x-rays 
from the radiation source. As the source 12 and the detector 16 are revolved about the examination 
region 14 via the rotating gantry 20, views are acquired over an angular range of preferably about 
360° or more. Optionally, a reduced scan of between about 180° and 360° is used. In another 
5 suitable embodiment, the detector 16 is replaced by a stationary detector ring mounted on a 
stationary gantry 22. Typically, a subject support 26 is linearly movable in an axial or z-direction 
that is generally transverse to a plane of rotation of the rotating gantry 20. 

Multiple-slice computed tomography projection data are suitably acquired by 
performing successive axial scans with the subject support 26 stationary during each axial scan and 

10 stepped linearly between axial scans. In this arrangement, the detector 16 can have either a single 
row of detector elements (that is, a one-dimensional detector) or a two-dimensional array of 
detector elements. Alternatively, helical computed tomography projection data is suitably acquired 
during continuous linear movement of the subject support 26 and simultaneous rotating of the 
gantry 20. This effects helical orbiting of the radiation source 12 relative to an imaging subject 

15 disposed on the subject support 26. A generally conical radiation beam and a two-dimensional 
radiation detector is preferably used for acquiring helical projection data. 

The outputs of detector elements of the radiation detector 16 are converted to 
acquired integrated attenuation projection values fid 0 that are stored in an acquired projection data 
memory 30. Each projection datum |xd 0 corresponds to a line integral of attenuation along a line 

2 0 from the radiation source 12 to a corresponding one of the detector elements of the detector 16. The 
projection data can be represented in a sinogram format in which each two-dimensional slice of the 
imaged region of interest is represented by a projection data array having coordinates of viewing 
angle (<J>) and line integral index (n). 

In the case of a fan-beam, cone-beam, or other original acquisition geometry 

2 5 having non-parallel rays, a rebinning processor 32 optionally rebins the projection data into parallel 

views. For a cone-beam geometry, such parallel-rebinned views typically include some uncorrected 
beam divergence in the cone-angle direction. Although parallel rebinning can improve 
computational efficiency, the artifact correction techniques described herein are typically most 
effective when performed in the original acquisition geometry, that is, without parallel rebinning. 

3 0 Hence, the parallel rebinning processor 32 is preferably omitted. 

For typical fan-beam and cone-beam geometries, the line integral index n suitably 
corresponds to a detector index indicating a detector element used to measure the projection of 
index n. It is contemplated, however, that the line integral index n may lack a direct 
correspondence with detector element number. Such a lack of direct correspondence can result, for 
3 5 example, from interpolation between rebinned proj ections. 
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A slice cycling processor 34 cycles through the sinograms corresponding to spatial 
slices and successively inputs each sinogram into an artifact-correcting reconstruction processor 40. 
The artifact-correcting reconstruction processor 40 performs a reconstruction in which artifacts 
introduced by high density regions, such as metal clips, high-density dental fillings, or the like, are 
5 substantially corrected, even for images that contain discontinuous high density regions or one or 
more medium density regions in addition to the one or more high density regions. The artifact- 
correcting reconstruction processor 40 makes use of a set of pixel density thresholds [HuJ 42 to 
classify pixels of the reconstructed image slice into at least a high density pixel class, a medium 
density pixel class, and a low density pixel class. 

10 With reference to FIGURE 2, in a preferred embodiment the set of density 

thresholds [HoJ 42 include five classes: a high density class designated "A" corresponding to high 
density regions such as dental fillings, metal implants, and the like; a medium density class 
designated "B" corresponding to bone or other medium density features; a low density class 
designated "C" corresponding principally to soft tissues; an air density class designated "E" 

15 corresponding to air pockets, ambient air surrounding the imaging subject, or the like; and a 
transition region designated "D tf corresponding to densities intermediate between the low density 
class "C M and the air density class "E M . A high density/medium density threshold H meto i defines a 
minimum density of the high density pixel class "A" and a maximum density of the medium 
density class "B". A medium density/low density threshold Hh C defines a minimum density of the 

20 medium density pixel class "B" and a maximum density of the low density class "C". A low 
density/transition density threshold Hj c defines a minimum density of the low density pixel class 
"C" and a maximum density of the transition density class ft D". A transition density/air density 
threshold Haitfmax) defines a minimum density of the transition density pixel class "D" and a 
maximum density of the air density class "E". 

2 5 Although in the preferred embodiment five density classes "A", "B", "C", "D", "E" 

are defined, it is also contemplated to omit the transitional density pixel class "D", in which case 
the low density/transition density threshold Hi c is omitted and the threshold H a u(max) also defines the 
minimum density of the low density pixel class "C". In this four-class embodiment, the threshold 
Hai^max) is optionally shifted to a higher density so as to divide the omitted transition pixel class "D" 

3 0 between the low density pixel class "C" and the air density pixel class "E". Still further, the air 

density pixel class "E" can also be omitted, leaving only the density pixel classes "A", "B", "C", 
with the transitional and air density pixel classes "D" and "E" subsumed into the low density pixel 
class "C". Moreover, the density classification system can include more than five density classes to 
provide improved density resolution. 
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With returning reference to FIGURE 1, for each input slice the artifact-correcting 
reconstruction processor 40 outputs a two-dimensional artifact-corrected reconstructed image. In 
multi-slice or helical compute tomography imaging, spatially successive artifact-corrected 
reconstructed image slices are accumulated in an image memory 46 to define a three-dimensional 

5 artifact-corrected reconstructed volume image. If, however, the acquired projection data is limited 
to a single slice of the region of interest, then the acquired projection data corresponding to the 
single slice is processed by the artifact-correcting reconstruction processor 40 and the image 
memory 46 stores a two-dimensional artifact-corrected reconstructed image. Optionally, projection 
data corresponding to one or more image slices are acquired over a selected time interval to provide 

10 a temporal series of artifact-corrected reconstructed image slices or image volumes representative 
of a temporal evolution of the region of interest. 

A video processor 50 processes some or all of the contents of the image memory 
46 to create a human-viewable image representation such as a three-dimensional rendering, a 
selected image slice, a maximum intensity projection, a CINE animation, or the like. The 

15 human-viewable image representation is displayed on a user interface 52, which is preferably a 
personal computer, a workstation, a laptop computer, or the like. Rather than or in addition to 
displaying an image representation, selected contents of image memory 46 can be printed on paper, 
stored in a non-volatile electronic or magnetic storage medium, transmitted over a local area 
network or the Internet, or otherwise processed. In a preferred embodiment, the user interface 52 

20 communicates with a computed tomography imaging scanner controller 54 to enable a radiologist 
or other operator to control the computed tomography imaging scanner 10 to construct an imaging 
session, modify an imaging session, execute an imaging session, monitor an imaging session, or 
otherwise operate the scanner 10. 

FIGURE 3 shows a block diagram of a preferred embodiment of the 

25 artifact-correcting reconstruction processor 40. The input is a projection data set P 0 60 
corresponding to a two-dimensional slice of the region of interest of the imaging subject. 
Preferably, the projection data set 60 includes projection data from a 360° revolution of the 
radiation source 12 about the examination region 14; however, reconstruction of a reduced 
projection data set providing, for example, about 180° angular coverage, is also contemplated. 

3 0 With reference to FIGURE 4, a diagrammatic sinogram representation 62 of the 

acquired projection data set P 0 60 is shown. The ordinate or y-coordinate of the sinogram 
representation 62 corresponds to view angle (|>, while the abscissa or x-coordinate of the sinogram 
representation 62 corresponds to line integral index n. The sinogram representation 62 includes an 
exemplary trajectory 64 of a high density region, and an exemplary trajectory 66 of a medium 

35 density region. In general, the slice can include an arbitrary number of high density regions and an 



- 6 - 



WO 2005/008586 



PCT/IB2004/002325 



arbitrary number of medium density regions, as well regions of other classes shown in FIGURE 2. 
The regions produces various corresponding sinogram trajectories that can cross one another. 

With returning reference to FIGURE 3, a filtered backprojection processor 70 
performs filtered backprojection of the projection data set P 0 60 to generate an uncorrected 
5 reconstructed image I 0 74. Although filtered backprojection is preferred, the processor 70 can 
alternatively implement substantially any type of image reconstruction algorithm that is compatible 
with the geometry of the projection data set P 0 60. 

With reference to FIGURE 5, as is known in the art, the presence of one or more 
high density regions in the slice typically causes the uncorrected reconstructed image Io 74 to 
10 include metal artifacts which generally manifest in the image as streaks extending away from the 
high density region. FIGURE 5 shows a diagrammatic representation of an exemplary uncorrected 
reconstructed image Io 74 which diagrammatically shows such streaks. The diagrammatic 
reconstructed images of FIGURE 5 are not intended to correspond to images of any particular 
anatomical structure, but rather diagrammatically illustrate features of the artifact-correcting 
1 5 reconstruction process that are typically observable in image space. 

With returning reference to FIGURE 3, pixels of the uncorrected reconstructed 
image Io 74 are classified by a pixel density classification processor 78 to generate a segmented or 
pixel density-classified image 80 in which pixel values are replaced by density classification 
index values corresponding to the high, medium, low, transition, and air density pixel classes of 
20 FIGURE 2. In a preferred embodiment, the pixel density classification processor 78 employs 
thresholding using the set of thresholds [HuJ 42 to classify each pixel of the uncorrected 
reconstructed image I 0 74 into an appropriate pixel density class. 

The pixel density-classified image Idass 80 is diagrammatically shown in FIGURE 
5, where regions of the image essentially consisting of pixels of a particular class are labeled by an 
25 appropriate class index selected from the density class indices "A", "B", "C", "D", "E n of FIGURE 
2. As indicated in FIGURE 5, the exemplary image includes a region of high density "A" which 
could be a metal implant and a region of medium density "B" which could be a region of bone, 
both contained within a region of low density "C" which could be soft tissue or the like. An annular 
region of transition density "D n surrounds the region of low density "C", and a region of air density 
30 "E" fills the periphery of the pixel density-classified image I class 80, corresponding for example to 
ambient air surrounding the imaging subject. 

Although the exemplary transition region in FIGURE 5 is an annular boundary 
region, the transition density class "D" is not limited to annular boundary regions, the thresholds 
H Ic and H a ii(max) defining the transition density class "D" are preferably selected so that the transition 
3 5 density class "D" includes air cavities and regions surrounding air cavities inside human subjects. 



- 7 - 



WO 2005/008586 



PCT/1B2004/002325 



For example, the transition density class ,f D" preferably includes air cavities in the head such as 
sinus cavities, the throat, nasal cavities, and the like. 

With returning reference to FIGURE 3, a density class averaging processor 84 
computes an average density value for pixels of the low density class "C". In preferred 
5 embodiments which include an air density class, the density class averaging processor 84 also 
optionally computes an average density value for pixels of the air density class W B". Averaging of 
the air density class is omitted in the embodiment that employs only three classes and omits the air 
density class altogether. 

A correction image processor 88 selectively replaces selected higher density pixels 

10 of the uncorrected reconstructed image I 0 74 with lower density values to produce a synthetic 
image 1^ 90. Specifically, pixels of the uncorrected reconstructed image Io 74 felling within the 
high density class "A" as indicated by the pixel density-classified image 1^ 80 are replaced by the 
average density value for pixels of the low density class "C". Similarly, pixels falling within the 
low density class "C" are replaced by the average density value for pixels of the low density class 

15 "C". Alternatively, the "A" and "C n class pixels can be replaced with a preselected value or 
spectrum of values which fall near the center of a nominal "C" density region. If an average density 
value for pixels of the air density class has been computed, pixels falling within the air density 
class "E" are suitably replaced by the average density value for pixels of the air density class "E". 
Pixels of the medium density class "B" are not replaced. If a transition density class "D" has been 

2 0 defined, then pixels of the transition density class "D" are also preferably not replaced. 

With reference to FIGURE 5, the synthetic image 1^ 90 is diagrammatically 
shown. It will be observed that the synthetic image I synt h 90 contains contrast principally due to the 
regions of medium density "B" and of transition density M D M . Regions of high density "A" are 
effectively removed from the synthetic image 1^ 90 by replacing pixels in these regions with the 
25 average density value for pixels of the low density class "C". Similarly, image contrast due to 
density variations in regions of low density "C" are substantially reduced by substitution of these 
pixels by the average density value for pixels of the low density class "C". 

With returning reference to FIGURE 3 and with additional reference to FIGURE 4, 
a forward projection processor 94 forward projects the synthetic image Isynth 90 to produce synthetic 

3 0 projection data P synt h 96. A diagrammatic sinogram representation 98 of the projection data set Psynth 

96 is shown in FIGURE 4. The sinogram representation 98 retains the exemplary trajectory 66 of 
the medium density region of the sinogram representation 62, since the pixels of medium density 
class "B" were not replaced. Similarly, trajectories due to regions of transition class "D" are 
retained in the synthetic projection data Psynth 96 since these pixels are retained in the synthetic 
3 5 image L^th 90. The sinogram representation 98 does not include the exemplary trajectory 64 of the 
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high density region of the sinogram representation 62, however, since the pixels of high density 
class "A" were replaced by the average density value for pixels of the low density class "C". 

In general, the sinogram representation 98 of the synthetic projection data set P^o, 
96 retains attenuation contributions to the line integrals due to regions of medium density class "B", 
5 even if those line integrals additionally have contributions from regions of high density class M A" in 
the acquired projection data set P Q 60. The forward projection processor 94 preferably projects the 
synthetic image Isynth 90 using a projection geometry corresponding to the geometry of the 
computed tomography imaging scanner 10. 

With continuing reference to FIGURE 3, the synthetic projection data set Pjynth 96 

10 is used for selectively replacing projection data of the acquired projection data set P 0 60 having 
attenuation contributions from high density regions. A high density region edge finding processor 
100 identifies edges of regions essentially consisting of pixels of the high density pixel class. 
Regions of high density can be identified, for example, by employing a moving analysis window 
and identifying high density regions as those having more than a selected number of pixels of the 

15 high density class within the analysis window. Once a high density region is identified, a suitable 
edge finding algorithm employs image differentiation or convolution to selectively enhance and 
identify edges. 

In an alternative approach, the edge finding processor- 100 performs binary 
thresholding on the uncorrected reconstructed image I 0 74 using the high density/medium density 

20 threshold Hmecal to produce a binary image having "l" binary values for pixels of the high density 
class and "0" binary values for other pixels. The binary "1" values are then filtered to remove 
outlying pixels of value "1" that have few or no neighboring pixels of value "1", and the remaining 
groups of pixels having value "1" define the high density regions. Edges are then identified in the 
filtered binary image as "©"-to-"!" binary value transitions, identified for example using an 

2 5 "exclusive-or" binary operator acting on neighboring pixel pairs. Those skilled in the art can readily 
employ other edge finding algorithms in constructing the edge finding processor 100. 

The edge finding processor 100 suitably outputs an edge-enhanced image lHD_edg« 
102 that identifies edges of high density regions. The edge-enhanced image iHD.cdge 102 is 
diagrammatically shown in FIGURE 5, where edges 106 of the region of the uncorrected 

30 reconstructed image Iq 74 consisting essentially of pixels of the high density class are indicated. In 
exemplary FIGURE 5, the edges 106 define a substantially circular edge; however, the high density 
regions can have substantially arbitrary shape with substantially arbitrary edges. Typically, the high 
density regions correspond to dental fillings, metal implants, or other compact discrete objects, and 
have well-defined edges defining simply-closed geometries. 
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With reference to FIGURE 3, a high density pixel replacement processor 110 
selectively replaces projections of the acquired projection data P 0 60 having contributions from 
high density regions with synthetic projection data from the synthetic projection data set P^u, 96 to 
produce a corrected projection data set Pi 112. In one approach, the high density pixel replacement 
processor 110 replaces projection data that intersect one or more high density regions identified by 
the edge-enhanced image ItiDedgo 102 with corresponding projection data from the synthetic 
projection data set P syn th 96. This approach, however, can lead to substantial attenuation 
discontinuities at transitions between the original acquired projection data P Q 60 and the 
replacement synthetic projection data set P synth 96. 

Hence, in a preferred embodiment, the high density pixel replacement processor 
110 interpolates between the acquired projection data and the synthetic projection data at the 
identified edges 106 of the high density image region. A suitable interpolating projection 
replacement formula replaces the acquired projection data of the projection data set P 0 60 with 
replacement projection data iid^x having values given by: 

Hd^ (n) = nd 8ymh (n) + + bf^O (1), 

where a=[fid 0 (ni)- jidsynthCnO], b=[fid 0 (n 2 )- ^d S ynth(n 2 )], index n is the line integral index, indices ni 
and n 2 are line integral indices of the edges 106 of the high density image region as shown in 
FIGURE 5, p.do indicates acquired projection data of the projection data set P 0 60, and nd^ 
indicates projection data of the synthetic projection data set P 8y „th 96. Projection data in the range 
n!<n<n 2 is replaced in accordance with Equation (1). Review of Equation (1) shows that 
Hdrepi(ni)= ^do(ni) and nd rep i(n 2 )~ |ido(n 2 ), providing a smooth transition at the edges 106. 

The interpolative replacement performed by the high density pixel replacement 
processor 110 operates on each view of the sinogram that includes one or more high density 
regions. That is, for each view specified by a given viewing angle (j), a given high density region 
has edges at line integral indices ni, n 2 which in general are different for different views of the 
same high density region. The interpolative replacement is performed for each view using the n X9 n 2 
edge values computed for that view. Moreover, there may be more than one non-contiguous high 
density region within a given view, with each such non-contiguous high density region having its 
own n x , n 2 edge index values. The interpolative replacement set forth in Equation (1) is repeated for 
each non-contiguous or non-overlapping high density region. 

With reference to FIGURE 4, a diagrammatic sinogram representation 114 of the 
corrected projection data set Pi 112 is shown. Advantageously, the sinogram representation 114 
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retains the trajectory 66 of the medium density region substantially intact, even where the trajectory 
66 of the medium density region crosses the location of the removed trajectory 64 of a high density 
region (indicated in the diagrammatic sinogram representation 114 by dashed lines representing the 
edges ni, n 2 of the replaced high density region in the views). There are no gaps of low density 
5 projections along the trajectory 66 of the medium density region because the synthetic projection 
data set P 8yn ih 96 retains contributions of medium density regions to the integrated attenuation 
projection values |id 0 while selectively removing contributions of high density regions to the 
integrated attenuation projection values jad^ Similarly, although not illustrated, a trajectory of a 
transition density region "D" remains intact even where such a trajectory is crossed by the high 

1 0 density trajectory 66. 

The corrected projection data set Pi 112 is input to the filtered backprojection 
processor 70 which performs filtered backprojection to generate an artifact-corrected reconstructed 
image 120. Preferably, a labeling processor 122 substitutes pixels defining a preselected label for 
pixels of the artifact-corrected reconstructed image 120 corresponding to high density regions that 

1 5 were substantially modified by the substitution of interpolated synthetic data. This labeling notifies 
a radiologist or other user who views a final artifact-corrected reconstructed image 124 that the 
labeled image regions are substantially modified by the artifact-correcting process. The pixels 
defining the preselected label can, for example, be pixels having high density values simulating the 
original high density region. In another approach, the pixels defining the preselected label define a 

20 selected pattern, such as a cross-hatch pattern, that is clearly artificial. The final artifact-corrected 
reconstructed image 124 is stored in the image memory 46 of FIGURE L 

The invention has been described with reference to the preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and understanding the 
preceding detailed description. It is intended that the invention be construed as including all such 

2 5 modifications and alterations insofar as they come within the scope of the appended claims or the 
equivalents thereof. 
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